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Abstract 

In this paper we construct global action-angle variables for the periodic 
Toda latticeQ 



1 Introduction 

Consider the Toda lattice with period N (N > 2), 

for n e Z, where the (real) coordinates {qn,Pn)nez satisfy {qn+N ,Pn+N) = 
{q-mPn) for any n € Z and the Hamiltonian HToda is given by 

N N 
n—1 n—1 

where a is a positive parameter, a > 0. For the standard Toda lattice, a ~ 1. 
The Toda lattice was introduced by Toda [TH] and studied extensively in the 
sequel. It is an important model for an integrable system of N particles in one 
space dimension with nearest neighbor interaction and belongs to the family of 
lattices introduce and numerically investigated by Fermi, Pasta, and Ulam in 
their seminal paper [S]. To prove the integrability of the Toda lattice, Flaschka 
introduced in [31 the (noncanonical) coordinates 

bn:=-Pn<^^, a„ := ae5(«"-«"+i) g M>o („ ^ g). 

These coordinates describe the motion of the Toda lattice relative to the center 
of mass. Note that the total momentum is conserved by the Toda flow, hence 
any trajectory of the center of mass is a straight line. 
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and the European Community through the FP6 Marie Curie RTN ENIGMA (MRTN-CT- 
2004-5652) 
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1 INTRODUCTION 



In these coordinates the Hamiltonian HToda takes the simple form 

N N 

2 



N N 



and the equations of motion are 



n=l 



n=l 



{n G Z). 



(2) 



Note that {bn+N, fln+Af) = (^Ti, On) for any n e Z, and n^=i = "^^^ Hence we 
can identify the sequences {hn)n& and (a„)„gz with the vectors {bn)i<n<N & 
and (a„)i<„<jv € lK>o- Our aim is to study the normal form of the system 
of equations ^ on the phase space 

M := X M^o- 

This system is Hamiltonian with respect to the nonstandard Poisson structure 
J = Jf,_a, defined at a point {b,a) = ((6„, a„)i<„<Ar by 



J = 



A 

-*A 



(3) 



where A is the 5-independent N x A^-matrix 

f ai ... 

-fli 02 
-02 03 



1 







V 



(4) 



-OAT-i OAT / 



The Poisson bracket corresponding to ([3]) is then given by 

{F,G}j{b,a) = {{VbF,VaF), J{VbG,VaG)U2. 

^ {WbF, A VaG)R« - {WaF, A* V,G)r« . 



(5) 



where F,G E C^(M) and where Vb and Va denote the gradients with respect 
to the A^-vectors b — (6i, . . . , ^at) and o = (oi, . . . , a^v), respectively. Therefore, 
equations can alternatively be written as bn — {bn,H}j^ d„ = {an,H}j 
{l<n< N). Further note that 



{bn,a7i}j — {bn+i,an}j — --^, 



(6) 



while {bn, Ofe} J = for any n, k with n ^ {k, k + 1}. 
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Since the matrix A defined by ^ has rank N — 1, the Poisson structure J 
is degenerate. It admits the two Casimir function^ 

Ci := J2 and C2 := I | | a„ ) (7) 

n— 1 

whose gradients Vfc^aCi = (VfcQ, VaCi) (i = 1,2), given by 

VbCi = V,Ci = 0, (8) 

VbC2 = 0, V„C2 - ^ f-, . . . , — V (9) 

N \ai QN J 

are linearly independent at each point (6, a) of M. 
The main result of this paper ist the following one: 

Theorem 1.1. The periodic Toda lattice admits globally defined action-angle 
variables. More precisely: 

(i) There exist real analytic functions {In)i<n<N-i on M which are pairwise 
in involution and which Poisson commute with the Toda Hamiltonian H 
and the two Casimir functions Ci, C2, i-c. for any 1 < m,n < N — 1, 

{IniJnjj^O OnM 

and 

{H, In} J = and {Cj, /„},/ = on M. 

(ii) For any 1 < n < N ^ 1 there exist a real analytic submanifold Dn of 
codimension 2 and a function ^„ : M \ Dn K, defined mod lix and real 
analytic when considered mod vr, so that on M, \ lJri=i^ ^"-r {^n)i<n<N-i 
and (In)i<n<N-i are conjugate variables. More precisely, for any 1 < 
TO, n < iV - 1, i 1,2 

{Im,0n}j = 5mn and {Ci,6'„}j = onM\Dn 



{^?m,^„}j = OnM\{D,n\jDn). 

Let M/3,a := {{b,a) e M^w . {CuC2) = {[3, a)} denote the level set of 
(Ci,C2) for a) € M x M>o. Note that (-/^Iat, al^v) € Mp^a where Ijv = 
(1, . . . , 1) g M^. As the gradients ^b,aCi and \/b,aC2 are linearly independent 
everywhere on M, the sets M/j^q are (real analytic) submanifolds of M of codi- 
mension two. Furthermore the Poisson structure J, restricted to M/j^q, becomes 
nondegenerate everywhere on M^g^o, and therefore induces a symplectic structure 
vi3,a on M/3,Q. In this way, we obtain a symplectic foliation of M with M/3,a 
being the symplectic leaves. 



smooth function C : M ^ R is a Casimir function for J if {C, ■},/ = 0. 
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Corollary 1.2. On each symplectic leaflVlp^a, the action variables {In)i<n<N-i 
are a maximal set of functionally independent integrals in involution of the pe- 
riodic Toda lattice. 



In subsequent work [5], we will use Theorem II. II to construct global Birkhoff 
coordinates for the periodic Toda lattice. More precisely, we introduce the model 
space y := x M x IR>o endowed with the degenerate Poisson structure 

Jo whose symplectic leaves are R^(^~^) x {/?} x {a} endowed with the standard 
Poisson structure, and prove the following theorem: 



Theorem 1.3. There exists a real analytic, canonical diffeomorphism 

n : (M, J) ^ (7, Jo) 

(5,a) 1-^ ((a;„,y„)i<„<Ar_i,Ci,C2) 

such that the coordinates {xn, yn)i<n<N-i, C'l , C2 are global Birkhoff coordinates 
for the periodic Toda lattice, i.e. {xmyn)i<n<N-i o,re canonical coordinates, 
Ci,C2 are the Casimirs and the transformed Toda Hamiltonian H — H o f2~^ 
is a function of the actions {In)i<n<N-i and Ci, C2 alone. 

In [TU] we used Theorem 11.31 to obtain a KAM theorem for Hamiltonian 
perturbations of the periodic Toda lattice. 

Related work: Theorem 11.11 and Theorem 11.31 improve on earlier work on the 
normal form of the periodic Toda lattice in [1] [2] . In particular, we construct 
global Birkhoff coordinates on all of M instead of a single symplectic leaf and 
show that techniques recently developed for treating the KdV equation (cf. 
[TT1[T2]) and the defocusing NLS equation (cf. [3 [16]) can also be applied for 
the Toda lattice. 

Outline of the paper: In section [2] we review the Lax pair of the periodic 
Toda lattice and collect some auxiliary results on the spectrum of the Jacobi 
matrix L(b,a) associated to an element (6, a) € M. In section [3] we study the 
action variables {In)i<n<N-i, and in section [4] we define the angle variables 
{()n)i<n<N-i on M \ U^^iDn usiug holomorphic differentials defined on the 
hyperelliptic Riemann surface associated to the spectrum of L(6, a). In sections 
[5] and [6] we establish formulas of the gradients of the actions and angles in terms 
of products of fundamental solutions and prove orthogonality relations between 
such products which are then used in section [7] to show that (/n)i<n<Af-i and 
{0n)i<n<N-i are canonical variables and to prove Theorem 11.11 and Corollary 
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2 Preliminaries 



It is well known (cf. e.g. [19]) that the system ([2]) admits a Lax pair formulation 
L — ^ — [B, L], where L = L^{b, a) is the periodic Jacobi matrix defined by 



L^ib,a) := 



(hi ai 
a2 ^3 



\±aN ■ ■ ■ 
and B the skew-symmetric matrix 

/ ai 

-fli 02 







B = 



-02 





few , 
-ajv \ 







(10) 



\aN ... 

Hence the flow of L = [-5, L] is isospectral. 

Proposition 2.1. For a solution (b{t),a{t)'^ of the periodic Toda lattice 
the eigenvalues {X^)i<j<N of L(b{t),a{t)^ are conserved quantities. 

Let us now collect a few results from 17 and 19J of the spectral theory of 
Jacobi matrices needed in the sequel. Denote by M"' the complexification of the 
phase space M, 

= {(6, a) e C^^ : Re > V 1 < j < N}. 

For (6, a) G M"' we consider for any complex number A the difference equation 

{Rb^avm = Xyik) {k e Z) (11) 

where y{-) = y{k)kez G and R^^a is the difference operator 

Rb,a ^ ak-iS-^ + bkS'' + OkS^ (12) 

with S"™ denoting the shift operator of order to e Z, i.e. 

(S""y)(fc) = y{k + to) for k e Z. 

Fundamental solutions: The two fundamental solutions yi(-,A) and y2{-,X) 
of (fTTj) are defined by the standard initial conditions j/i(0, A) = 1, yi(l, A) = 
and 2/2(0, A) = 0, 2/2(1, A) = 1. They satisfy the Wronskian identity 



W{n) yi(n, X)y2{n + l,X)-yi{n + 1, A)y2(ri, A) 



ON 



(13) 
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Note that for n = N one gets 

W{N) ^ 1. (14) 

For each k G N, yi{k, \,b,a), i = 1,2, is a polynomial in A of degree at most 
k — 1 and depends real analytically on (b, a) (see T71). In particular, one easily 
verifies that y2{N + 1,A, 6, a) is a polynomial in A of degree N with leading 
coefficient a^^ . 

Wronskian: More generally, one defines for any two sequences (w(?i))„ez and 
(w(n))„gz the Wronskian sequence (W^(^^))„gz = ^f')("))nez by 

W{n) :— v{n)w{n + 1) — v{n + l)w(n). 

Let us recall the following properties of the Wronskian, which can be easily 
verified. 

Lemma 2.2. (i) If y and z are solutions of kll]) for A = Ai and A — A2, 

respectively, then W — W{y, z) satisfies for any /c G Z 

akW{k) = ak-iW{k - 1) + (A2 - Xi)y{k)z{k). (15) 



(a) If y{-, A) is a 1 -parameter- family of solutions of ill]) which is continuously 
differentiable with respect to the parameter A and y{k, A) :— ■^y{k, A), 
then W — W{y, y) satisfies for any k G li 

akW{k)^ak-iW{k-l) + y{k,\f. (16) 

Discriminant: We denote by A(A) = A(A, 6, a) the discriminant of (jlip . 
defined by 

A(A) :=j/i(iV,A)+2/2(iV + l,A). (17) 

In the sequel, we will often write Aa for A(A). Note that y2{N + 1,A) is 
a polynomial in A of degree N with leading term a~^A^, whereas yi{N,X) 
is a polynomial in A of degree less than N, hence A(A, 6, a) is a polynomial 
in A of degree N with leading term a~^A^, and it depends real analytically 
on (6, a) (see e.g. [in])- According to Floquet's Theorem (see e.g. JB,), for 
A G C given, pT|) admits a periodic or antiperiodic solution of period N if the 
discriminant A(A) satisfies A(A) = +2 or A(A) = —2, respectively. (These 
solutions correspond to eigenvectors of i+ or i^, respectively, with defined 
by pop .) It turns out to be more convenient to combine these two cases by 



(a) Af = 4 

(b) AT = 5 

Figure 1: Examples of the discriminant A (A) 



considering the periodic Jacobi matrix Q — Q{h, a) of size 2N defined by 



Q = 



bi 


ai 

















ai 


b2 





















a-N- 


-1 bN 


ajM 


















UN 


bi 


ai 

















ai 


b2 






a-N 

















1 bN 



Then the spectrum of the matrix Q is the union of the spectra of the matrices 
and L~ and therefore the zero set of the polynomial — 4. The function 
A^ — 4 is a polynomial in A of degree 2N and admits a product representation 

2N 

Al-A^a-'^^lliX^X,). (18) 

The factor a^^^ in (fT8|) comes from the above mentioned fact that the leading 
term of A(A) is a-^A^. 

For any {b, a) € M, the matrix Q is symmetric and hence the eigenvalues 
i^j)i<j<2N of Q are real. When listed in increasing order and with their alge- 
braic multiplicities, they satisfy the following relations (cf. [17]) 

Al < A2 < A3 < A4 < As < . . . X2N-2 < X2N-I < A2Ar. 

As explained above, the Xj are periodic or antiperiodic eigenvalues of L and 
thus eigenvalues of L+ or according to whether A(A) = 2 or A(A) — —2. 
One has (cf. [H]) 

A(Ai) = (-1)^-2, A(A2„) = A(A2„+i) = (-l)"+^-2, A(A2iv) = 2. (19) 
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2 PRELIMINARIES 



Since is a polynomial of degree N, A\ = A(A) = ^A(A) is a polynomial 
of degree A'^ — 1, and it admits a product representation of the form 



N-l 



Na-"" lliX-Xk). (20) 



k=l 

The zeroes (A„)i<„<Ar_i of A^ satisfy A2n < A„ < X2n+i for any I < n < 
N — I. The intervals (A2n, A2n+i) are referred to as the n-th spectral gap and 
7n '■— A2n+i — A2n as the n-th gap length. Note that |A(A)| > 2 on the spectral 
gaps. We say that the rt-th gap is open if 7„ > and collapsed otherwise. The 
set of elements (6, a) E M for which the rt-th gap is collapsed is denoted by Z3„, 

A. :={(6,a)eM:7„ = 0}. (21) 

By writing the condition 7„ = as 7^ = and exploiting the fact that 7^ 
(unlike 7„) is a real analytic function on M, it can be shown as in 12J that Dn 
is a real analytic submanifold of M of codimension 2. 

Isolating neighborhoods: Let (b, a) G M be given. The strict inequality 
A2n-i < A2n guarantees the existence of a family of mutually disjoint open 
subsets {Un)i<n<N-i of C so that for any 1 < n < — 1, J7„ is a neighborhood 
of the closed interval [A2n, A2ri+i]. Such a family of neighborhoods is referred to 
as a family of isolating neighborhoods (for (5, a)). 

In the case where {b,a) G M"', we list the eigenvalues {Xj)i<j<2N in lexico- 
graphic ordering 

Ai A2 ^ A3 ^ . . . ^ A2Ar. 

We then extend the gap lenghts 7„ to all of M"' by 

7n := A2„+i - A2„ (1 < n < TV - 1) 

and define 

i?^={(&,a)e W:7„ = 0}. (23) 

In the sequel, we will omit the superscript and always write Z?„ for D^. 

Similarly, we do this for the zeroes (A„)i<„<Ar_i of A^. The A^'s and A^'s 
no longer depend continuously on (6, a) G M^. However, if we choose a small 
enough complex neighborhood W of M in M^, then for any {b, a) G W the closed 
intervals G„ CC (1 <n<A^— 1) defined by 

G„ {(1 - t)X2n + tX2n+i : < t < 1} (24) 

are pairwise disjoint, and hence, as in the real case, there exists a family of 
isolating neighborhoods {Un)i<n<N~i- 
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The lexicographic ordering a ^ b for complex numbers a and b is defined by 

Re a < Re b 

a^b -.-^^ { or (22) 
Re a = Re b and Im a < Im b. 
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Lemma 2.3. There exists a neighborhood "W o/ M in M*^ such that for any 
(6, a) € W, there are neighborhoods Un of Gn in C (1 < n < N ^ 1) which are 
pairwise disjoint. 

Remark 2.4. In the sequel, we will have to shrink the complex neighborhood 
W several times, but continue to denote it by the same letter. 

Contours r„; For any (6, a) € W and any 1 < n < — 1, we denote by r„ 
a circuit in [/„ around G„ with counterclockwise orientation. 

Isospectral set: For (&, a) G M, the set Iso(6, a) of all elements (6', a') £ M 
so that Q{b', a') has the same spectrum as Q{b, a) is described with the help of 
the Dirichlet eigenvalues /ii < /12 < • • ■ < Mtv-i of (jlip defined by 



yi{N + l,^i„) = 0. 



(25) 



They coincide with the eigenvalues of the (N— 1) x (N— l)-matrix L2 = L2{b, a) 
given by 

/ 62 a2 ... \ 



0-2 





Vo 



On 





on-i 
^1 bN I 



In the sequel, we will also refer to fii, . . . , fj-N-i as the Dirichlet eigenvalues 
of L{b,a). Evaluating the Wronskian identity (fT3|) at A = /x„ one sees that 
/i„ lies in the closure of the n-th spectral gap. More precisely, substituting 
yi(iV + 1, = in the identity (fT3|) with A = /i„ yields 

(26) 



yi(iV,^„)y2(iV + l,/i„) = 1. 
Hence the value of the discriminant at fi„ is given by 



MPn) = y2{N + I, + 



1 



y2(^ + l,Atn) 



(27) 



and |A(/x„)| > 2. By Lemma [2.61 below, given the point {b,a) with bi = . . . = 
bN — P and ai ~ ... ~ apf ~ ol, one has A2n = A2„+i and hence /i„ = A2n 
for any l<n<A^— 1. It then follows from a straightforward deformation 
argument that A2n ^ fJ-n l£ A2n+i everywhere in the real space M. 

Conversely, according to van Moerbeke [17], given any (real) Jacobi matrix 
Q with spectrum Ai < A2 < A3 < A4 < A5 < . . . A2JV-2 < -^2^-1 < A2Ar 
and any sequence (/i„)i<„<Ar-i with A2„ < fJ-n < >^2n+i for n 1, . . . , iV - 1, 
there are exactly 2^ A^-periodic Jacobi matrices Q with spectrum (A„)i<„<2Af 
and Dirichlet spectrum (/iri)i<n<Af-i, where r is the number of n's with A2n < 

/J-n < A2n+1- 

In the case where (6, a) € M.^, we continue to define the Dirichlet eigenvalues 
{lJ-n)i<n<N-i by (j25p . and we list them in lexicographic ordering fii -< /i2 -< 
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... ^ fiN-i- Then the ^i's no longer depend continuously on (6, a) G M"'. 
However, if we choose the complex neighborhood W of M in "M"" of Lemma [231 
small enough, then for any (&, a) G W and 1 < n < iV — 1, /i„ is contained in 
the neighborhood J7„ of G„ (but not necessarily in G„ itself). 

Riemann surface Sb^o-' Denote by Sb^a the Riemann surface obtained as the 
compactification of the afRne curve Gb,a defined by 

{{X,z) eC^ : = A^{X,b,a)~A}. (28) 

Note that Cf,^a and Eb^a are spectral invariants. (Strictly speaking, "Sb.a is a 
Riemann surface only if the spectrum of Q{b, a) is simple - see e.g. Appendix A 
in [H] for details in this case. If the spectrum of Q{b,a) is not simple, I](6, a) 
becomes a Riemann surface after doubling the multiple eigenvalues - see e.g. 
section 2 of [13].) 

Dirichlet divisors: To the Dirichlet eigenvalue /i„ (1 < n < N — 1) we 
associate the point on the surface Sb.a, 

:= (;u„, ^A2^ -4] with 4^~l ■■= 2/i(iV,A*n) - y2(A^+ l,Ai„), (29) 
where we used that, in view of (pS]) and the Wronskian idcntiy ([T^ . 

^L-4 = (yi(^.Mn)-2/2(iv + i,/^„))'. 

Standard root: The standard root or s-root for short, \/l — A^, is defined for 
AeC\[-l,l]by 

\/l - A2 := ^A Vl - A-2. (30) 

More generally, we define ior \ £ C \ {ta + {1 — t)b\0 < t < 1} the s-root of a 
radicand of the form (6 — A) (A — a) with a < b, a =/= b hy 

</ib-X)iX-a) ^ i/l-w^, (31) 

where j := b ~ a, t := and w := 

Canonical sheet and canonical root: For (6, a) S M the canonical sheet of 
Eb_a is given by the set of points (A, ^ — 4) in Cb,a, where the c-root A^ — 4 

is defined on C \ U^=o('^2n, A2n+i) (where Ao := —00 and A2Ar+i :— 00) and 
determined by the sign condition 



-i^A2-4>0 for A2Ar_i < A < A2JV. (32) 

As a consequence one has for any 1 < 71 < 

sign ^A2_^„-4= (-1)^+"-! for A2„ < A < A2„+i. (33) 

The definition of the canonical sheet and the c-root can be extended to the 
neighborhood W of M in M"' of Lemma 12.31 
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The s-root and the c-root will be used together in the following way: By the 
product representations pO)) and p8|) of Ax and A| — 4, respectively, one sees 
that for any {b, a) in W \ with 1 < n < iV — 1, 



iV(A-A„) 



V^I^ \/(A2„+i-A)(A~A2„) 



X„(A) VAer„ (34) 



where 



( — 1)^+"^-^ A — Am 

^"^^^ " V(A-Ai)(A2^^ V(A-A2™+i)(A-A2™)' ^^^^ 

Note that the principal branches of the square roots in ([35|) are well defined for 
A near Gn and that the function Xn is analytic and nonvanishing on Un- In 
addition, for (b, a) real, Xn is nonnegative on the interval (A2n, A2n+i)- 

Abelian differentials: Let (b,a) € M and 1 < n < iV — 1. Then there exists a 
unique polynomial ipnW of degree at most A^— 2 such that for any 1 < fc < A^— 1 

Here, for any l<fc<A^— 1, Cfe denotes the lift of the contour Tk to the 
canonical sheet of Sfc.a- For any k ^ n with A2A; =/= ^2k+i, it follows from (j36p 
that 

1 /■^^'' + ^ V'n(A) 



TT 



VAT 



: dA = 0. (37) 



Hence in every gap (A2fc, A2fc+i) with k ^ n the polynomial tpn has a zero which 
we denote by cr^. If A2fc = A2/C+1 then it follows from p6|l and Cauchy's theorem 
that CT^ = A2fc — A2fc+i. As '0n(A) is a polynomial of degree at most N — 2, one 
has 

V'n(A)=M„ H {\-al), (38) 



l<fc<Af- 



where M„ = M„(6, a) ^ 0. 

In a straightforward way one can prove that there exists a neighborhood W 
of M in Mc, so that for any (6, a) € W and any 1 < n < N —1, there is a unique 
polynomial V'n(A) of degree at most N —2 satisfying ([55]) for any 1 < k < N — 1 
as well as the product representation (|38[) . and so that the zeroes are analytic 
functions on W. 

Lemma 2.5. Let 1 < n < N — 1 be fixed. Then the zeroes {Xk)i<k<N-i of 
A(A) and {<7^)i<k<N-i,k^n of ipnW satisfy the estimates 

Xk-rk= 0{jI), (39) 
- Tfe = Oi^l). (40) 

near any given point {b, a) G W, where — ■h{X2k+i + A2A:). 
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Proof. To verify write - 4 in the form 

- 4 = (A - A2„)(A2„+i - AK(A) (41) 



where p„ is a polynomial which does not vanish for A G C/„ . Then (|39p follows 
by differentiating (|4T|l with respect to A at A„. 

Fix 1 < A:, n < A^— 1 with k ^ n. In a first step we prove that a'^—Tk — 0{'-fk) 
near any given point (b,a) E W. If -fk ~ 0, then cr^ = r^, and (|40p is clearly 
fulfilled. Hence we assume in the sequel that 7^ 0. By the product formulas 
p8|) and p8)) for '0n('^) ^nd A| — 4, respectively, we obtain, for A near Gk, 



^A^ v/(A2fe+i -A)(A-A2fe)' 
where 

(A - <) V(A - Ai)(A2jv - A) „|j,, V(A2m+i - A)(A2^ - A) ' 

with a" :— Tn and M^^{b,a) ^ 0. The function (J! is analytic and nonvanishing 
in Uk- Substituting (|42|) into ([36]) one gets 

' =Cr(A)dA = 0. (44) 



27^^rfc v/(A2fe+i - A)(A - A2fe) 

We now drop the superscript n for the remainder of this proof and write (^k as 
a(A) = a + (Cfc(A) - a) with a- := Cfc(n-) 7^ 0. Note that 

1 /■ A-CTfe 

a\ = Tk — (Jk 



7r;, v/(A2fe+i - A)(A - A2fc) 
and hence ((44|) becomes 

r u W (A-^fc)(a(A)-a) 

2^ -'r, V(A2fc+i - A)(A - A2fe) 

To estimate the integral on the right hand side of (|45p . note that 

/(A) 



27r 7rfc </(A2A;+i - A)(A - A2fc) 



■.d\ 



< max |/(A)| (46) 



for an arbitrary function / analytic on Uk- We want to apply (|46|) for /(A) = 
(A - afc)(Cfc(A) - a). Note that for A e Gfc, 

ia(A)-ai = ia(A)-a(Tfc)i < Af|7fei, 

where M = suplJ^<j.<^_;^{|Cfe(A)| : A e Gfc}. Hence (gHl) leads to 
kfc - TfellCfcl = sup |A - (Tfcl 0(7fc). 
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Dividing by |^fc| ^ 0, we get 



\cTk - Tk\ = sup |A - CTfcl 0{-fk) 



(47) 



and in particular \(Jk — Tk\ — 0{'-fk)- 

In a second step, we now improve the estimate (H71) . Note that 



sup I A - o-fcl < IcTfe -Tk\+ sup I A - Tfcl = 0(7fe). 



(48) 



Combining (IT71) and (I48p . we obtain the clained estimate (HDl) . 



□ 



For later use, we compute the spectra of Q{b,a) and £2(6, a) in the spe- 
cial case (6, a) = (/JIat, alAr) with /3 e M and a > 0. Here Ijv denotes the 
vector (1, . . . , 1) G M^. These points arc the equilibrium points (of the restric- 
tions) of the Toda Hamiltonian vector field (to the symplectic leaves M/j^q). We 
compute the spectrum {Xj)i<j<2N of the matrix Q{(3lN,alN) and the Dirich- 
let eigenvalues {^ik)i<k<N-i of L — L(/31jv, alAr) together with a normalized 
eigenvector gi = (5Ki))i<j<Ar of A^i> i-^- Lgi = mgi, gi{l) = 0, and a vector 
hi = (^i(i))i<j<7v '^liich is the normalized solution of Ly — pny orthogonal to 
gi satisfying W{hugi){N) > 0. 

Lemma 2.6. The spectrum {Xj)i<j<2N ofQ{[31iq, OilN) and the Dirichlet eigen- 
values (/i/)i<;<7v-i of L{(31n ,ol1n) are given by 



In particular, all spectral gaps of Q{P\n ,Oi\N) are collapsed. For any 1 <l < 
N — 1. the vectors gi and hi defined by 



X21 — X21+1 — f^i 
X2N 



(3 -2a, 

/3-2acos— (l</<iV-l), 
13 + 2a. 




(49) 



(50) 



satisfy Ly = fJ-iy and the normalization conditions 



N N 



Y.gi{jf^Y.^i{jf^l, gi{0)>0, <?;(!) = 0; 



W{hi,gi){N)>0, {hi,gi)^.=0. 
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2 PRELIMINARIES 



Remark 2.7. For {h, a) — {PIn, aljv) the fundamental solutions yi and y2 are 
given by 

mUA)--'^^^^^^^^ ijez) (5i) 

smp 

mU,X) = '^ Uez) (52) 

sm p 

where n < p < 2Tr is determined by cos p = . 

Proof. For any A G R, the difference equation pT|) for (/SIat, aljv) reads 

{Rp.ay)(k) f3y{k) + ay{k - 1) + ay{k + 1) = Ay(fc) (53) 
and can be written as 

y{k-l) + y{k + l) = ^^y{k). (54) 
a 

Since we are looking for periodic solutions of (|54p . we make the ansatz y{k) = 
Q±ipk_ rp]_^jg igads to the characteristic equation 

2cosp=e'P + e-'P = ^— ^. 

a 

For the solution to be 2Ar-periodic, it is required that p E jjZ. To put the 
eigenvalues in ascending order, introduce pi = {1 + j^)tt with < I < N . Then 
for any 1 < j < 2N, there exists < I < N such that 



Xj = P + 2a cos pi = (3 — 2a cos — 



In 
N' 



Note that for Z = 0, Ai = /3— 2a is an eigenvalue of Q{f3lN, al m) with eigenvector 
y{k) = e^'^'' — (— l)*^- Similarly, for I — N, = /3 + 2a is an eigenvalue with 
eigenvector y{k) = 1. For the eigenvalue X21 = f3 — 2a cos ^ {1 < I < N — 1), 

y±{k) = e^'P^'^ 

are two linearly independent eigenvectors. As there are 2N eigenvalues allto- 
gether, X21 is double for any 1 < I < N — 1, and Ai and X2N are both simple. It 
follows that all A^ — 1 gaps are collapsed and hence pi = X21 for all 1 < I < A^ — 1. 

Turning to the computation of and hk, one easily verifies that for any real 
number A ^ ±2a + /3, the fundamental solution yi{-, A) of ([51)1 with yi{0, A) = 1 
and ?/i(l, A) = is given by 

sin (p(j — 1)) 
sm p 

where tt < p < 27r is determined by cos p — , thus proving (|5ip . In the same 
way, one verifies ([52]) . For X — pi — 13 — 2a cos ^ we then get 

sin {piij 1)) = sin ((1 + l)^(j - 1)) = (-1)^+1 sin ^^^il^. 
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In particular, sin — 1)) = for j = I and j = N + 1. As 



and these two sums add up to N, one sees that 

N 



2^sin ^ = (55) 

yielding the claimed formula (j49|) for g;. 

By the same argument one shows that hi given by (— l)-' y^cos ^''^^^''^ (i-e. 
the right side of ([50]) ) satisfies Rp,ahi = fiihi and the normahzation condition 
"^jLi — 1- Using standard trigonometric identities one verifies that 

N 

and W{hi, gi){N) can be computed to be 

hi{N)giiN + l)-hi{N + l)gi{N) = -hiiN + l)<?i(iV) = ~hi{l)giiO) > 0. 

Hence hi is indeed the eigenvector with the required normalization, i.e. hi = hi, 
thus proving ([50)1 . □ 

3 Action Variables 

In the next two sections, we define the candidates for action-angle variables 
on the phase space M and investigate some of their properties. In this sec- 
tion we introduce globally defined action variables (/ri)i<n<A'-i proposed by 
Flaschka-McLaughlin [4J . 

Definition 3.1. Let (6, a) G M. For 1 < n < iV - 1, 

= dX (56) 



27r7r„ ^Af 

where A\ — -j^^x is the X-derivative of the discriminant Ax = A(A, 6, a) and 
the contour r„ and the canonical root ^ are given as in section\^ 

Remark 3.2. The contours r„ can he chosen locally independently of {b,a). 
In view of the fact that A\ is a spectral invariant of L(h,a) the actions In are 
entirely determined by the spectrum of L{b,a). In particular, {In)i<n<N-i are 
constants of motion, since by Proposition \2. 1\ the Toda flow is isospectral. 
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3 ACTION VARIABLES 



Remark 3.3. The variables {In)i<n<N-i can also be represented as integrals 
on the Riemann surface 'Sb,a ■ If Cn denotes the lift of r„ to the canonical sheet 
o/E^^o, 156]) becomes 

In = :^f A r^- dX (l<n<iV-l). (57) 



From the definition (|56p . the foUowing resuh can be deduced: 

Proposition 3.4. On the real space M, each function /„ is real, nonnegative, 
and it vanishes if "in =0. 



Proof. Since 



it follows that 



/„ = i-/^^(A-A„)^^rfA. (58) 



By shrinking the contour of integration to the real interval, we get 



In = - / (-l)^+"-i(A - A„) ^,^,; dX 



by taking into account the definition ([32]) of the c-root. Since sign(A — A„)Aa — 
(—1)^+"^^ on [A2ri, A2n+i] \ {A„}, the integrand is real and nonnegative, hence 
/„ is real and nonnegative on M, as claimed. 

If 7„ = 0, then A2n = A2n+i. Hence A„ = A2n = A2„+i = r„ and 



A - A„ = i \/ (A2„+i - A)(A - A2„). 

Therefore the integrand in (]56p is holomorphic in the interior of the contour r„, 
and by Cauchy's theorem the integral in (|55]) vanishes. □ 

The action variables (/n)i<ri<jv-i can be extended in a straightforward way 
to a complex neig hborhood W of M in M"^. 

Theorem 3.5. There exists a complex neighborhood W of ?vi in such that 
for all 1 < n < N — 1, the functions In defined by i56\) extend analytically to 
W, /„ : W ^ C. 

Proof. Let W denote a neighborhood of M in M"' of Lemma 12.31 and define for 
any 1 < n < — 1 the functions /„ on W by the formula (|56| . Let {b, a) e W 
be given. Then there exists a neighborhood Wf,.a of (b, a) in W so that the 
integration contours r„ in (]56p can be chosen to be the same for any element 
in Wfc^a and Aa/\/A^ — 4 is analytic on B^{Tn) x Wb.a, where i?e(r„) := {A € 
C| dist (A, r„) < e} is the e- neighborhood of r„ with e sufficiently small. This 
shows that /„ is analytic on W. □ 
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Proposition 3.6. There exists a complex neighborhood W of JVl in M"' such 
that for any 1 < n < N ^ 1, the quotient In/jn extends analytically from 
JA \ Dn to all o/W and has strictly positive real part on W. As a consequence, 
£,n = ^/'^In/ln '^^ analytic and nonvanishing function on W, where ^ is 
the principal branch of the square root on C\ (— oo, 0]. 



Proof. Let W be the complex neighborhood of Theorem 13.51 Substituting (|34p 
into ([58|l leads to the following identity on W \ Dn 

where x„ is given by ([35)) . Upon the substitution X{() = Tn + with t„ = 



^{^2n + ^2n+i) and Sn = i o^i^ then obtains 



2(A„-r„) 



7^ 4Wr;yW^^"(^" + f^)'^^' ^''^ 

where T'^ is the puUback of r„ under the substitution A = A(^), i.e. a circuit in 
C around [—1, 1]. By A„ — = 0(7^), and hence (5„ ^ as 7„ ^ 0. We 
conclude that 

2/„ N f . . CdC , .N t^dt N 



lim„ ^ -r X«(t„)-7=== = X„(t„)— / = -rXniTn) 



072 47r7r;/^ ^^/W^ ^"^"'2^7-1^1^ 



N 

x«(r„)-y_^ 



4 



By defining 2^ by -jXn ('''«) on W n -D„, it follows that is a continuous 
function on all of W. This extended function is analytic on W \ _D„ as is its 
restriction to WnZ3„. By Theorem A. 6 in [15] it then follows that "^f- is analytic 
on all of W. 

By Lemma 13.71 below, the quotient Inlln *^an be bounded away from zero 
on M, ^ > ij^^j^iixj- shrinking W, if necessary, it then follows that for 
any 1 < n < TV — 1, the real part of In/in positive and never vanishes on W. 
Hence the principal branch of the square root of 2J„/7^ is well defined on W 
and ^„ has the claimed properties. □ 

To show that \J^^ is well defined on W, we used in the proof of Proposition 
13.61 the following auxiliary result, which we prove in Appendix [Al 

Lemma 3.7. For any {b, a) e M and any 1 < n < N — 1, 

ll<M^2N -\l)In- (60) 

From the definition (|56p. Proposition [331 and the estimate ([60| one obtains 
Corollary 3.8. For any (6, a) G M and any 1 < n < N — 1, 

In = if and only if 7„ — 0. 
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4 ANGLE VARIABLES 



Actually, Lemma [3. 71 can be improved. We finish this section with an a priori 
estimate of the gap lengths 7„ in terms of the action variables and the value of 
the Casimir C2 alone, which will be shown in Appendix B. 

Theorem 3.9. For any (6, a) e M^_q with /3 e M, a > arbitrary, 

N-l /N-l \ /N~l \ ^ 



n— 1 \n— 1 



4 Angle Variables 

In this section, we define and study the angle coordinates {9n)i<n<N-i- Each 
9n is defined mod 27r on W \ £)„, where W is a complex neighborhood of M in 
as in Lemma 12.31 and Dn is given by ([23|) . 

Definition 4.1. For any 1 < n < N — 1. the function On is defined for {b, a) € 
M\Dnby 

N-l 

dn ■■^Vn+ E ^fc mod2Tr, (62) 

n^k=l 

where for k ^ n, 

PI ^ r^^dA, ^„ = r-^W dA [mod2.l (63) 

a77.d where for 1 < k < N ^i*f, is the Diriehlet divisor defined in I129\) . and \2k 
is identified with the ramification point (A2fc,0) on the Riemann surface Sfc.a. 
The integration paths on S^^a in i6S\) are required to be admissible in the sense 
that their image under the projection n : Eb,a ^ C on the first component stays 
inside the isolating neighborhoods Uk- 

Note that, in view of the normalization conditions (|36|) of ip„, the above 
restriction of the paths of integration in (|63p implies that rjn and hence On are 
well-defined mod 27r. 

Theorem 4.2. Let W be the complex neighborhood of M, in M*^ introduced in 
Lemma \2.3l Then for any 1 < n < N ^1, the function On '■ W\D„ — )■ C {modn) 
is analytic. 

Remark 4.3. As the lexicographic ordering of the eigenvalues ofQ{b,a) is not 
continuous on W, it follows that rjn and hence On are only continuous mod n on 
W. 

Proof of Theorem \4.2\ To see that dn : 'W\ D„ C (mod tt) is analytic, define 
for any 1 < fc < — 1 the set 



Ek := {{b,a) £ : Hk{b,a) € {A2fe(6, a), A2fc+i(6, a)}}. 
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Below, we show that for any 1 < A; < TV — 1 with k ^ n, (3]!: is analytic on 
W \ (Dfe UEk), that its restrictions to Dfc n W and EkHW are weakly analyticQ, 
and that it is continuous on W. Together with the fact that EkHW and Dk H W 
are analytic subvarieties of W it then follows that is analytic on W - see 
Theorem A. 6 in [12]. Similar results can be shown for /3" = ry„ (mod tt) on 
W \ Dn, and one concludes that 9n (mod tt) is analytic on W \ _D„. 

To prove that fc 7^ n, is analytic on W \ {Dk U Ek), note that since X2k 
is a simple eigenvalue on W \ Dk, it is analytic there. Furthermore, /i^ is an 
analytic function on the (sufficiently small) neig hborhood W of M in M^. On 
W \ {Dk U Ek) we can use the substitution A = A2fe + z to get 

A,, VaT^ Jo V^^/W) 

where D{z) — ^ (Mk+z)-4: analytic near z — and D{0) 7^ 0. Note that 
D{z) does not vanish for z on an admissible integration path not going through 
A2fe+i. Such a path exists since (6, a) is in the complement of Ek- Furthermore 
'ipn{X2k + z) and D{z) are analytic in z near such a path and depend analytically 
on (6, a) e W\ (DfcUE'fe). Combining these arguments shows that /3|f is analytic 
on W\ {DkUEk). 

For k n with A2fc 7^ A2/C+1 one has 

"^"(^^ -.dX = 0. (64) 



A2 



"'fe ~ ^2A: if A2fe = A2A;+i One sees that continues to hold for {b, a) e 
EkHW with A2fe = A2fc+i and we have PklEk^w = 0. To prove the analyticity 
of PklokHW consider the representation ^ of -^J^^- For {b,a) e Dk n W, 
one has 

A2fe = A2fc+1 — Tk — (Tk, 

which implies that the factor , ^ '^k jj-^ (j42|) equals ±i. Hence we 

V (A2fc + i— A)(A— A2fc) 

can write 

'■'^'^ ^„(A) rf^^ 



PI - / dX^±z C^{X)dX. 

h,, v/A^-4 ^ 

As /ifc is analytic on W, it then follows that Pkln^nw is analytic. To see that 
is continuous on W, one separately shows that is continuous at points in 
W\{DkUEk), EkHWXDk, DknW\Ek and DkOEkDW, where for the proof 
of the continuity of at points in Dk O Ek D W we use (|^^ and the estimate 
- Tfc = 0{-fl) of Lemma [231 

By (j64p . r/„ vanishes mod tt on En H W \ _D„. Arguing in a similar way as 
for one then concludes that ry„ (mod tt) is analytic on W \ _D„. □ 



*Let E and _F be complex Banach spaces, and let U G E be open. The map / : J7 — > F is 
weakly analytic on U, if for each u £ U , h £ E and L G F*, the function 2 1— > Lf{u + 2/1) is 
analytic in some neighborhood of the origin in C. 
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5 GRADIENTS 



5 Gradients 

In this section we establish formulas of the gradients In, On {1 < n < N ~ 1) 
on M in terms of products of the fundamental solutions yi and ?/2- 
Consider the discriminant for a fixed value of A as a function on M, 

Axib,a)=yi{N)+y2{N+l). 

Then A\ is a real analytic function on M. To obtain a formula for the gradients 
of yi{N) and y2{N + 1) with respect to b, differentiate Rb,ayi = ^yi with respect 
to b in the direction v G to get 

{Rb^a - X){Wby^,v){k) = ~vky^{k). (65) 

Differentiating Rb,ayi = Mji with respect to a in the direction u G leads to 

(i?b,a - \){Vayi,u)(k) = ~Uk~iyi{k - 1) - Ukyi{k + 1). (66) 

Taking the sum of ([55]) and (pS)) yields 

{Rb,a - A) ((Vf,2/^, v) + (Va2/z, u)) (fc) = -{Rv,uyz){k) (67) 

which we can rewrite as 

(i?b,a- A)(V6,ay^, («,«)) (fc) = -{R,,uy^){k), (68) 

where (•, •) in ((68|) now denotes the standard scalar product in R^^, whereas in 
(|65|) . ([66]) . and ([67|) it is the one in M.^ . The inhomogeneous Jacobi difference 
equation ([68|) for the sequence {Vb^ayi, {v,u)){k) can be integrated using the 
discrete analogue of the method of the variation of constants used for inhomo- 
geneous differential equations. As (Vb^a2/i, (w, u))(0) — (Vb^ayi, (w, u)) (1) = 0, 
one obtains in this way for m > 1 



(Vb.a yi, {v, u)){m) = - ( y^^'^^ yi{k){Ry^uyi){k) 



yi{m) 



Y,y2ik){Rv,uyim]. (69) 



In the sequel, we will use (j69l) to derive various formulas for the gradi- 
ents. The common feature among these formulas is that they involve products 
between the fundamental solutions yi and y2 of Whereas the gradients 

with respect to 6 = (6i, . . . , 6jv) involve products computed by componentwise 
multiplication, the gradients with respect to a = (oi, . . . , a^r) involve products 
obtained by multiplying shifted components, reflecting the fact that the bj are 
the diagonal elements of the symmetric matrix Lib, a), whereas the Oj are the 
off-diagonal elements of L{b, a). 
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To simplify notation for the formulas in this section, we define for sequences 
{v{j)jei.), {w{j)j<^7,) Q C the iV- vectors 

vw:= (w(i)w(j))i<^.<^ , (70) 

vSw := {v{j)w{j + l))i<j<Ar , (71) 



where S denotes the shift operator of order 1. Combining ([70]) and ([71]) . we 
define the 2A^-vector 

v-sw := {v ■ w,v ■ Sw + w ■ Sv). (72) 
In case v = w we also use the shorter notation 

:= v-sv. (73) 
Written componentwise, v-sw is the 2A^-vector 



(w-sw)(j> 



v{j)w{j) (1 < i < N) 

v{j-N)w{j-N+l) + v{j-N + l)w{j-N) {N<j <2N) 



Proposition 5.1. For any {b,a) € M, the gradient \/b,a^\ — (V^A^, VqA^) 
is given by 

-ajvVbAA - y2{N)yi-yi - yi{N+l)y2-y2 + {y2{N+l) - yi{N)) yi-y2{74:) 
-aNVaAx = 2y2{N)yi ■ Syi - 2yi{N + l)y2 ■ Sy2 

+ {y2{N + 1) - yi{N)) [yi ■ Sy2 + • %i) (75) 

or in the notation introduced above 
Vb.aAx = - — (y2(iV)yf ~yi{N + l)yi + (y2(iV + 1) - y,{N)) 2/1-^^2) • 

ON 

(76) 

The gradients V^Aa and VaA^ admit the representations {I < m < N) 

— — 2/2(iV,A,5™6,^"a), (77) 
^ = - f — y2(A^+l, A, ^"6, 5™a) + ^y2(A^-l, A, S^'+\ S'''+^a)\ .(78) 

00,m \am Ctm+l / 

Proof. The claimed formula ([7S)) follows from the dcfintion of A^ and formula 
Indeed, evaluate ([55)1 for i = 1 and to = to get 



N 

Wb,ayi, («,-«)) (iV) = -^^^V 2/1 (fc) (wfc_i?/i(fc - 1) + vkyiik) + ukyi{k + 1)) 
On ^-^ 

k=l 

+^^^y2y2{k) {uk-iyi{k - 1) + Vkyiik) + ukyi{k + 1)) . (79) 
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5 GRADIENTS 



In order to identify these two sums with {yf, (w, u)) and (yi-s2/2, ('i', respec- 
tively, note that 



N N 

Uk-iyi{k)yi{k - 1) = X] ^kyi{k)yi{k + 1) + matTi 
fc=i fc=i 

where 

:=j/i(0)yi(l)-j/i(iV)2/i(iV + l). 

For the second sum in Ij79p . we get an expression of the same type witli a similar 
correction term 

T2 := yi(0)y2(l) - yiiN)y2iN + 1). 

Taking into account the initial conditions of the fundamental solutions and the 
Wronskian identity p^ . one sees that y2{N)Ti — yi{N)T2 vanishes. Hence we 
have the formula 

(Vb,a2/i, iv,u))iN) = -— {y2{N){vl, {v,u)) - yi{N){y^-sy2, («,«))) • (80) 

Oat 

Similarly, evaluating formula for i = 2 and m — N + 1 leads to 

{Vb,ay2Av.u)){N + l) = -— {y2{N+l){yi-sy2.{v,u))^yi{N+l){yl,{v,u))). 

(81) 

Here we used that the value of the right side of ([M)) does not change when we 
omit the term for fc = TO = A'^-|-lin both sums. 

It remains to prove the two formulas ([77)1 and ((78|) . We first note that 



y2{n, X, S"^b, S^^a) — (y2{n + m, \,b, a)yi{m, X,b, a) 
aN V 

-yi{n + Tn,X,b,a)y2im,X,b,afj, (82) 



since both sides of ([5^ are solutions of Rs"^b,S"^ay — Xy (for fixed m £ Z) 
with the same initial conditions at n = 0, 1. For n — 1 this follows from the 
Wronskian identity (jl3|) . Similarly, one shows that 



yi{N + m.,X) = yi{N,X)yi{m,X)+yi{N + l,X)y2{m,X), (83) 
y2{N + m,X) = y2iN,X)yiim,X)+y2{N + l,X)y2im,X). (84) 

for any (6, a) G M. Hence, suppressing the variable A, we get 
y2{N, S^b, S'^a) ^ ({y^{N)y,{m) + y2{N + l)y2(m)) y,{m) 

UN ^ 

-{yi{N)yi{m) + yi(iV + 1)?;2(to)) y2{m)^ 
= —(y2{N)yi{mf + {y2{N+l)~ y^{N))yi{m)y2{ni) 
~yi{N + l)y2imf^ 
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By ^ this leads to 

and formula ((77|) is established. To prove ([78|) . we first conclude from (|82|) that 
-^y2(Af - 1, S'^+^fe, S-^+^a) = y2{N + m, 6, a)yi(m + 1, 6, a) 

+ TO, 6, a)y2(™ + 1, a) (85) 

and 

— y2(^ + 1, ^"a) = V2{N + TO + 1, 6, a)yi(TO, 6, a) 

—yi{N + m+l,b,a)y2{m,b,a). (86) 

Now expand the right hand sides of ([85]) and ((86| according to ((83)) and ([84)) . 
By dZS]), the sum of ([85]) and §Q is -ajvf^, thus proving ([78]). □ 

As a next step, we compute the gradients of the Dirichlet and periodic eigen- 
values. In the following lemma, we consider the fundamental solution ?/i(-,/i) 
as an A^-vector yi(i, /i)i<j<JV- Let \\yi{fi)\\'^ = Tlf=i 2/1 (JiM)^, and denote by ' 
the derivative with respect to A. 

Lemma 5.2. If fi is a Dirichlet eigenvalue of L{b,a), then 

aNyliN,^l)yl{N + I, ^^) = ||2/i(m)II' > 0. (87) 

In particular, yi{N + 1, /i) =/= 0, which implies that all Dirichlet eigenvalues are 
simple. 

Proof. This follows from adding up the relations (fTB)) . □ 

As the Dirichlet eigenvalues {^in)l<n<N-l of L{b,a) coincide with the roots 
of yi (A^-l- 1, /i) and these roots are simple, they are real analytic on M. Similarly, 
the eigenvalues Ai and A2Ar are real analytic on M, whereas for any 1 < n < 
A^ — 1, A2ti and A2,i+i are real analytic on M\ £>„. Note that for {b, a) G M\Dn 
and i £ {2n, 2n + 1}, we have A^^ 7^ as A.; is a simple eigenvalue. 

Proposition 5.3. For any 1 < n < N ~ 1, the gradients of the periodic eigen- 
values Xi (i = 2n, 2n + I) on M. \ Z)„ and of the Dirichlet eigenvalues fin on M 
are given by 

Vfe.aAi = ^= ^ y2 Vb,al^n = 9n, (88) 

where we denote by fi the eigenvector of L{b,a) associated to Xi, normalized by 

N 

^/,(j)2 = l and (/,(!), /,(2)) e (M>oxM)U({0}xM>o), 
i=i 

and where = i9nij))i<j<N is the fundamental solution j/i(-,/i„) normalized 
so that Yl!j=i9n{i f = 1- 
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5 GRADIENTS 



Proof. We first show the second formula in ([55)1 . Differentiating yi(7V+l, = 
with respect to {b, a), one obtains 

yi(A^ + l,/i„) 

Here we used that yi{N + 1, /i„) 7^ by Lemma [5.21 To compute the gradient 
'^b,ayi{N + 1, A)|a=;^„, we evaluate ((69|) for i = 1 and m = N + 1. In view of 
j/i(A^ + 1, fin) — and taking into account (|26|) . one then gets 

aNyi{N, fin)yi{N + l,/i„) 

The claimed formula 'Vb,afJ-n — dn t^ow follows from Lemma 15.21 By differenti- 
ating Aa; — ±2 with respect to (6, a), one obtains Vb,a\ — —^b,a^\\\=\i/^\i 
in a similar fashion. To see that ^b,a^i — fi, differentiate Rb,afi — Kfi with 
respect to (6, a) in the direction {v,u) £ R^^, 

RbA^b,af^Av,u)){k) + {R,^uf^){k)^{\7b,aX^Av,u))Mk)+M^b,af^Av,u)){k)^ 

where (•, •) denotes the standard scalar product in M?^ . Take the scalar product 
(in M.^) of the above equation with fi. Now use that 

{^b,afi{v,u),Rb,afi)K« (Vf,,a/i («, w) , 

(/j,/<)r" = 1, and 

{Rv^ufi, fi)M'^ = {fi ,(w,u))r2n, 

to conclude that S/b^aK — ft holds. □ 

To compute the Poisson brackets involving angle variables we need to estab- 
lish some additional auxiliary results. Recall from section [3] that for 1 < k,n < 
N ~ 1 with k n and (6, a) G M, is given by 

,;.r^.x, (91) 

whereas 

A" := Vn = r 4== d\ (mod 2n). (92) 



By Theorem 14.21 the functions (3^ with k ^ n are real analytic on M, whereas 
/3^\ when considered mod tt, is real analytic on M \ Z)„. 

Proposition 5.4. Let 1 < k < N — 1 and (6, a) G M. //7fc > and \2k = l^k, 
then for any 1 < n < N ~ 1, 

^b,aPk = T gk-slT-k, 
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where hk denotes the solution of Rb_ay ~ tJ-kU orthogonal to g^, i.e. 

N 

J29kij)hk{j)^0, 
satisfying the normalization condition W{hk, gk){N) = 1. 

Proof. We use a limiting procedure first introduced in [16j for the nonlinear 
Schrodinger equation and subsequently used for the KdV equation in [Tl] , [12] . 
We approximate {b,a) G M with X2k{b,a) = ^k{b,a) < A2fc+i(6, a) by (b',a') G 
Iso(6, a) satisfying A2fc(6, a) < ^k{b\a') < \2k+i{b,a). For such {b',a'), using 
the substitution A = A2fc + z in the integral of (|9ip . we obtain 

Pkib', a') = r 44^ . r'-' M^,,, (93) 



where D{z) = D{X2k, z) := (A^(A2fe + z) — 4)/z. Taking the gradient, undoing 
the substitution, and recalling the definition ([29]) of the starred square root then 
leads to 

Vfc,a/?^ = f"^^'^ (Vb,a/ifc - V,,aA2fc) + E{b\ a'), (94) 
/A2 ^4 

with the remainder term E{b' , a') given by 

Eib',a')^ r%,, fMX.k + .)\dz 



^D{X2k,z) J ^A 

As the gradient in the latter integral is a bounded function in z near z = 0, 
locally uniformly in {b',a'), it follows by the dominated convergence theorem 
that lim(f,/^a')-»(b,a) E{b' ,a') = 0. 

The gradient Vb,aPk depends continuously on (6, a) E M, hence we can 
conclude by (|94p that it can be written as 



Vb,al3'^= lim f"^^'^ (Vfc'.g'Mfc- Vfc',a'A2fc). (95) 

(6',a')^(b,a) ^•/A2 - 4 



The gradient of both sides of the Wronskian identity (fl3|) . 

2/1 (TV, A)y2(iV + 1, A) - yi(iV + 1, X)y2{N, A) = 1, 

leads to 

yiiN + l)Vb,ay2{N) + y2{N)Vb,ayiiN + 1) (96) 
^y2{N+l)Wb^aA + {yi{N)-y2{N + l))Vb,ay2{N+l), 

The A-derivative can then be computed to be 

yi{N + l)y2{N) + yi{N + l)y2{N) 

=y2{N + 1)A + {yi{N) - y2{N + \))y2{N + 1). (97) 
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5 GRADIENTS 



Using (IMl), dSZD, and yi{N + 1, fik) = 0, formula ((Ml) for Vt,a^J.k leads to 
2/2(iV + l)Vb,aA + (yi(iV) - y2{N + l))Vb,ay2iN + 1) 



Further by 



y2{N + 1)A + {yi{N) - y2{N + l))y2{N + 1) 

Vb,aA 



A 



A=A2 



(98) 



(99) 



Now substitute ^ and ^ into {S/h'^a'fJ^k - Vh',a'A2fc) and use that by (pS)) . 



^/a2~4 = (yi(iV) - 2;2(iV + 1))^,. 

We claim that 

Vb^aAifc-Vb^aAafc _ 2/2(iV+l)Vf,,ayi(A^)-2/i(A^)V6,ay2(iV+l) 



lim 

(!>',"') 
^(6, a) 



Ayi(7V + l)y2(A^) 



(100) 



Indeed, to obtain (|100p after the above mentioned substitutions, we split the 
fraction , ,1 . .a'l^-k — '^b',a'^2k) into two parts which are treated sepa- 



V f'fc 



rately. In the first part we collect all terms in — (S/b' a'l^k ~ Vf,' a'^2k) 
which contain [yi{N) — y2{N + 1))!^^ in the nominator, 

-AU,, • yb',a'y2{N + i)U, + yb',a'A\x,, ■ miN + i)U, 



/(a', 6') 



AU,, • iy2{N + 1)A + (yi(7V) - y2iN + l))y2iN + 
Using again ((97)) we then get 

UN + l)Vb,ayiiN) - yiiN)Vb,ay2iN + 1) 



lim I(b , a ) 

(b' ,a')^(b,a) 

The second term is then given by 



A2/i(iV + l)y2(A^) 



AU,, • (yi(iV) - y2(iV + 1))|,,, • (z/i(iV + 1)2;2(^))U/ 

Note that the nominator of II{b',a') is of the order 0(/ifc — A2fe)- In view of 
([29ll . we have 

(yiiN) - y2{N + 1))U, = 0(V/^,-A2fe) 

whereas the other terms in the denominator of II{b' , a') are bounded away from 
zero. Indeed, A2fc being a simple eigenvalue for (5, a) means A 1^21- 7^ for (fe', a') 
near (6, a). Further, use a version of ((97l) in the case A2fe — Hk to conclude that 

yi(A^ + i)y2(iV) = y2(A^ + i)AA2,. 

Hence yi{N + l)y2{N)\^^ 7^ for {b',a') near (6, a) and II{b',a') vanishes in 



the limit of 



^2fc- 
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Substituting ((BD]) and into ^UU^, we obtain 
miN + l)Vb,a yiiN) - yi{N)Vb,a y2{N + 1) 



Ayi(JV + l)y2(iV) 

!/l(JV+l) ''^ 



Mfc Oat A 



1 



with 2/0 = ^"(j^+n j/i - ?/2- Hence 



Vb,a/3fe = — -yrsyo- 

Since is invariant under the translation b i— > b+t{l, .... 1), the scalar product 
(Vb^a/3fe , (1, 0))ji2N vanishes. Hence 



It means that yi and yo are orthogonal to each other. Finally we introduce 
hk '■= \\yi\\ya and verify that 

W{hk,gk) = W fbibo, - W{yo,yi) = ('M^±|lyi - y^, 

V \m\\J \yiiN + i) 

By ([131), it then follows that 

W{hk,gk) = -Wiy2,yi) =W{yi,y2). 

Hence by pi)) . 



M^(/ife,5fc)W = W^(yi,Z/2)W = l. 
This completes the proof of Proposition l5.4l □ 



6 Orthogonality relations 

In Propositions 15.11 15. 3[ and 15.41 we have expressed the gradients of A\, fin, 
and, on a subset of M, of /3^f in terms of products of fundamental solutions of 
the difference equation pTjl . In this section we establish orthogonality relations 
between such products - see [2j for similar computations. Recall that in ([7^ we 
have introduced for arbitrary sequences (uj)jgz, (w^jOjez the 2A^-vector v-sw. 

Lemma 6.1. For any (6, a) G M, let vi, wi and V2, W2 be pairs of solutions of 
ill]) for arbitrarily given real numbers fi and A, respectively. Then 

^^^'^^ {V1-SWI,J{V2-SW2))=V + B, (101) 

aiON 

where 

V := {Wi- SW2)\^ + {SWi-W2)\^ (102) 
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6 ORTHOGONALITY RELATIONS 



with Wi and W2 denoting the Wronskians Wi := W{vi,W2), W2 '■— W{wi,V2), 
and where B is given by 

B-=^^^ {{vi ■ {V2-SW2){2N) ~ {V2 ■ ws) If {v^-swi){2N)) . (103) 

Proof. We prove poip by a straightforward calculation, using the recurrence 
property of the Wronskian sequences Wi and W2- By the definition ^ of 
J we can write 

2 {V1-SWI,J{V2-SW2)) ^ Ei+ Bi, 

where 

N 

Ei:=Y,ak [{vi-sWi){k){v2-sW2){N + k) - {vi-swi){k + l){v2-sW2){N + k) 

k=l 

-{vi-sWi){N + k){v2-sW2)(k) + {vi-sWi){N + k)(v2-sW2){k + 1)] 

and 

B, :^ajv {{vrsWi)iN + 1) - (vi •su'i)(l)) («2-sW2)(2A^) 

+aN {{V2-SW2){1) - {V2-SW2){N + 1)) {vi-sWi)(2N). 

Let us first consider Ei. Calculating the products Vj-^Wj according to ([7^ . we 
obtain, after regrouping, 

N 

Ei^J2 ak[{v2{k)wi{k) +V2{k + l)wi{k + l))Wi{k) 
fc=i 

+ {vi{k)w2{k) +vi{k + l)w2{k + l))W2{k)] 
+B2 

with 

B2 a^iviiN + l)wi{N + 1) - (wi-sWi)(iV + 1)) (wa -s W2)(2iV) 
+aN{{v2-sW2){N + l) ~ V2{N+l)w2{N + l)) {vi-sWi){2N). 

Multiply El by (A — /i) and use the recurrence relation (fTSl) to express (A — 
/i)w2(fc)wi(fc), (A - /x)u2(fc + l)wi{k + 1), (A - ii)vi{k)w2{k), and (A - fj,)vi{k + 
l)w2{k + 1) in terms of the Wronskians Wi and W2 to get 

N 

(A - = 5]K«fc+i(W^i(fc)M^2(fc + 1) + W^i(fc + l)M^2(fc)) 
fc=i 

-afc_iafe(VFi(fc - l)VK2(fc) + Wiik)W2ik - 1))] 
+ {X-fi)B2. 

The sum on the right hand side of the latter identity is a telescoping sum 
and equals the term aiaj^V with V defined in (I102|) . In a straightforward way 
one sees that ^^(Si + B2) equals the expression B defined by (|103[) . hence 
formula (|10ip is established. □ 
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Corollary 6.2. For any A,/i G C, 

{Aa,A^}j = 0. (104) 
Proof. By the formula ([75]) for the gradient of Aa, 

{Aa,A^}j = {Vb,aAx,JVb,aAf,) 

is a linear combination of terms of the form (wi -s wi, J(u2 -s ^2)) for pairs of 
fundamental solutions vi, wi and V2, W2 of (fTTj) for /j, and A, respectively. In 
view of ([77)1 and ([75)) . V;,Aa and VqA^ are both A^-periodic. In the case A 7^ ^ 
we use Lemma I01 and note that the boundary terms (|102p and (|103l) in Lemma 
16.11 vanish, hence {A^, A^}j — 0. In the case A = /i the identity p04p follows 
from the skew-symmetry of {•, -jj. □ 

Corollary 6.3. Let I < k < 2N and A G C. On the open subset of M where 
Afe is a simple eigenvalue of Q{b,a) one has 

{Afe,AA}j = 0. 

Proof. Using formula ([551) f^'' '^b,a^k, we conclude from CoroUarv 16.21 that 



□ 



{Xk, Aa}j = -^{A^, Aa}jU=a. = 0. 



Corollary 6.4. Let be the n-th Dirichlet eigenvalue of L{b, a) and A 7^ 
a reaZ number. Then 

(A-^„)(yf(/i„), Jy?(A)) = faA.^i^^^±ii^y (105) 

a ^/ 2/ ^ , n\ 2 ?/i(iV + 1, A)j/2(iV + 1, A) 
(A-^„)(?/i(^„), Jyi(A)-sy2(A)) = -, ^2 ^^'^^> 



(A-..0«(.,.).J»|(A)).<,5,(^(^|f±i^J -Ij (107) 

Proof. The three stated identities follow from Lemma 16. 1[ using that yi {N + 
l,/i„) = 0, j/i(2,/i„) — — flAr/ai, and, by the Wronskian identity yi(iV,/x„)- 
yzCA^ + l,/^™) = 1. □ 

Corollary 6.5. Let fin be the n-th Dirichlet eigenvalue of L{b,a) and A =/= fin 
a real number. Then 



2 



;/i(iV + l,A) , ... 

{/Zn, AaIj = . , 1 7—; • 108) 

yi(A' + l,^„) A-/i„ 
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6 ORTHOGONALITY RELATIONS 



Proof. By ^U^, combined with (gS]), we get 

{/x„, A,}^ = y'^^"'^\^'^\ {y^,{f,^),JV,^aAx). (109) 

aNyi(N + I, fin) 

Substituting the formula ([75)) for JVfc.aA^ we obtain 

{i/, i^ln ) , JVfc,a A a) =- — 2/2 (iV, A) {yf {^in ) , Jyf (A)) 
aN 

- — (2/2 (iV + 1 , A) - yi (TV, A)) (j/f (Ai„) , Jyi (A) -s 2/2 (A)) 
aN 

+—yiiN + 1, A)(2/f (/x„), J2/f (A)). (110) 

Oat 

To evaluate the right side of (jllOp . we apply Corollary 16.41 and get 

^(2/?(Mn), JVmAa)= 7T^^- ^ ( - 2/2(^, A)2/i(7V + 1, A)2 

aiOAT ai2/2(iV + l,Mn)^ 

-(2/2(iV + 1, A) - 2/1 (TV, A))2/i(TV + 1, A)2/2(TV + 1, A) 

yi{N + l,X) 



+2/i(TV+l,A)2/2(TV + l,A)2) 



ai 



Using the Wronskian identity (|14p . the sum of the terms in the square bracket 
of the latter expression simplifies, and one obtains 

A - M» , 2( . . . yi{N + l,X) y,{N + l,X) 

aiaN ai2/2(iV + 1,M«) ai 

-(2/i(TV,A^„)2-l), (111) 



2/i(iV+l,A), ,2 



ai 

where for the latter equality we again used p4)). Substituting piip into (|109p . 
we get 

A - M» r A 1 2^2(A^+ l,M«)2/i(^ + l.A) 2 1^ 

aiO-N aiaMyi(N + 1, ^n) 

yiiN + l,\) 



l/S? - 4 
aiaAr2/i(A^+ l,Ai„) V ^" 

where we used that, by the definition of the starred square root (P^)) . 



= yi{N, ^„) - y^iN + 1, Ai„) = 2/2(^ + 1, ^in){yl{N, ^lnf - i). 

This proves (fTM)) . □ 
Proposition 6.6. For ani/ AeM, l<rt<iV-l, anrf (6, a) e M \ D„, 

{0„,AA}j=l/'n(A). (112) 
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Proof. Recall that 9,, = Y^k^i (mod 27r) with (3^ given by To 
compute {/3^,Aa}j, we first consider the case where (6, a) ^ V}k=i ^k and 
A2/C < Mfc < -^2^+1 for any 1 < fc < — 1. Then A2fe and /LtJ are smooth near 
(6, a) and, by Leibniz's rule, we get 



{jUfe, Aa} J , {A2fc, Aa}j 



>VA^, - 4 ^Ai^^ - 4 

By Corollary 16.31 {A2A:,Aa}j = 0. Moreover, as the gradient Vb,a ^^^^^4 
is orthogonal to a Iso (5, a) and JVb.a^x € Th.a Iso (6, a) it follows that the 
Poisson bracket { ^^2^\ ' Aa}j vanishes for any /i in the isolating neighborhood 
Un of Gn. Hence 

{/3^AA},/ = 4^k^{Mfc,AA},. 



A?, 



By (fTOS]) . we then obtain 



fr[ yi{N + l,fik) A - ^fe 

where for the latter equality we used that ^^=1 i)i(7v+i''^fc) a^^I '^^ 

are both polynomials in A of degree at most N ~ 2 which agree at the N — 1 

points {tik)i<k<N^i- 

In the general case, where (&, a) G M \ £>„ and the Dirichlet eigenvalues are 
arbitrary, A2fc ^ fJ-k ^ A2fc+i for any 1 < fc < iV — 1, the claimed result follows 
from the case treated above by continuity. □ 

Proposition 6.7. Let 1 < n,m, k,l < N — 1 and let (6, a) e M with X2i{b, a) = 
^ii{b, a) for i = k,l. Then 

{(3Lf3r}j = 0. 



Proof. In view of Proposition l5.4[ this amounts to showing that the scalar prod- 
uct {{gk'shk), Jigcshi)) vanishes. For k — I, this follows from the skew-symmetry 
of the Poisson bracket, hence we can assume k ^ I. We apply Lemma |6. II with 
vi := gk, wi := hk, W2 hi and W2 := gi, which implies that Wi = W{gk,gi) 
and W2 = W{hk, hi). Since 5^(1), gi{l), gk{N + 1) and gi{N + 1) all vanish, 
we conclude that VFi(iV) = 1^1(0) = and iSWi){N) = {SWi){0) = 0, hence 
the expressions V and E, defined in (|102p and (|103[) . vanish. This proves the 
claim. □ 
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7 CANONICAL RELATIONS 



7 Canonical relations 

In this section we complete the proof of Theorem 11.11 and Corollary 11.21 In 
particular we show that the variables (/n)i<n<w-i, (^ri)i<n<A^-i satisfy the 
canonical relations stated in Theorem ll.il 

Using the results of the preceding sections, we can now compute the Poisson 
brackets among the action and angle variables introduced in section [31 

Theorem 7.1. The action-angle variables {In)i<n<N-i o,nd {0„)i<n<N-i sat- 
isfy the following canonical relations for I < n,m < N — I: 

(i) on M, 

{/n,/m}j = 0; (113) 

(a) on M \ Dn, 

{On,Im}.J = —{Im,()n\j ~ —^nrn- (114) 



Proof. Recall that arcosh (t) = [t^ — 1) 2 . Hence for any (6, a) e M 



dt 



2tt ,lr d\ 



and therefore 

27r dX - 4 



Integrating by parts we get 

Vb.aAA 



27r Jr^ v/A^ 



dX. (115) 



As {Aa, A^},7 = for aU A, e C by CoroUarvlOl it follows that {/„, /„},/ = 
on M for any 1 < n, m < A'^ — 1. 

To prove (|114p . use pisp and then Proposition 16. 61 to get 

{0„,A;,}j^^_ 1 f M^) ^x = -6^ 



{(^nJ..h = -^J^^^^—dX = --J^ 



nm 1 



v/Aj^ 2^ Jr™ ^Al^ 

by the normalizing condition (j36p oi ipn- CH 

To prove that the angles (6'n)i<ri<Af-i pairwise Poisson commute we need 
the following lemma. We denote hy K = K{b, a) the index set of the open gaps, 
i.e. 

K{b, a) = {1 < n < iV - 1 : 7„(6, a) > 0}. 
Lemma 7.2. At every point (6, a) in JA, the set of vectors 

(l) ((Vb,a/„)nei<-:Vb,aCi,V6,aC2) 

and 
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are both linearly independent. 

Proof. The claimed statements follow from the orthogonality relations stated 
in Theorem 17.11 Let (&, a) G M and suppose that for some real coefficients 
(j'n)neK C R and si, S2 G K we have 

rn'^b^aln + SlWh^Cl + S2VC2 = 0. 

For any m £ K, take the scalar product of this identity with J\7b.a(^m- Using 
that {Imdm}j = 5nm and that Ci and C2 are Casimir functions of {•,•},/ one 
obtains 

Thus = for all m £ K, and it follows that si'^b,aCi + S2Vb,aC2 = 0. By 
([8|) and Vb,aCi and are linearly independent, hence si = S2 = 0. 

This shows (i). The proof of (i) also shows that (ii) holds. □ 

Theorem 7.3. In addition to the canonical relations stated in Theorem \7. 1\ the 

angle variables {On)i<n<N-i satisfy for any 1 < n,m < N — 1 on M\ {DnU D„i) 

{en,Om}j = 0. (116) 

Proof. Let 1 < n,m < N — 1. By continuity, it suffices to prove the identity 
(|116p for (&, a) e M \ {{jf^i^ Di). Let (6, a) be an arbitrary element in M \ 

{{jf^i^ Di) ■ Recall that Iso(5, a) denotes the set of all elements {b',a') in M 
with spec(Qh',a') = spec(Qfc,a), 

Iso (6, a) = {{b\ a') e M : A(-, 6', a') = A(-, b, a)}. 

Then Iso(&, a) is a torus contained in M \ ( -^0 ' ^^'^ eigenvalues of 

Q{b, a) are simple, its dimension is iV — 1. By Lemma [7?^ at any point (6', a') G 
Iso(6, a), the vectors {JVb',a'Ik)i<k<N-i are linearly independent. Using the 
formula (|115p for the gradient of Ik, one sees that, by CoroUarv 16. 2i for any 
/leR, l<fc<7V-l, 

Hence for any {b',a') e Iso(6, a), 

(J Vfc',Q'/fc)i<fc<Ar_i e Tb'^a' Iso (6, a), 

and therefore these vectors form a basis of Tb'.a' Iso (6, a). 
To prove the identity (|116p . we apply the Jacobi identity 

{F, {G, H}j}j + {G, {H, F}j}j + {H, {F, G}j}j = 
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to the functions Ik, On and 0„j. Since by Theorem 17. 1[ {Ik,On}j = Skn, we 
obtain 



{/fc, {On, 0m}j}j = on M \ IJ A for any 1 < fc < - 1. 



. 1=1 



It then follows by the above considerations that 'Vb',a'{dn,dm}j is orthogonal 
to Tb'^a' Iso (b, a) for all {b' , a') € Iso(&, a), i.e. {6'„, 6m} j is constant on Iso(&, a), 

{0n, OrajAb', a') - {0„, 0„J j(6, a) V (&', a') S Iso (6, a). 

By [17], Theorem 2.1, there exists a unique element (&', a') g Iso(&, a) satisfying 
IJLk{h',a') = \2k{h,a) for all 1 < fc < iV - 1. The claimed identity (fTTe]) then 
follows from Proposition 16. 71 □ 



Proof of Theorem \l.l\ By Theorem l3.5l and Theorem 14. 2[ the action and angle 
variables introduced in Definitions 13.11 and 14. li respectively, have the claimed 
analyticity properties. The canonical relations among these variables have been 
verified in Theorem 17.11 and Theorem 17.31 and the relations {Ci,/„}j = (on 
M) and {Ci,6n}j = (on M \ I?„) follow from the fact that Ci and C2 are 
Casimir functions. It remains to show that the actions Poisson commute with 
the Toda Hamiltonian. To this end note that that the Hamiltonian H can be 
written as 

N N N 

H-\T.^l + T.-l = \'^ (Lib, af) = \ Y^iX^f 

n— 1 n— 1 j=l 

where {^^)i<j<N are the N eigenvalues of L{b,a). Recall that on the dense 
open subset M \ U^SiDk of M, the A^'s (1 < i < N) are simple eigenvalues 
and hence real analytic. It then follows by (|115p that for any 1 < n < — 1, 



where for the latter identity we used Corollarv l6.3l Hence for any 1 < n < N~l, 

{i/,/„}j = onM\U^~^^Dk. 

By continuity it then follows that {H,In}j — everywhere on M. □ 

Proof of Corollary ] LSI Since for any /? G M and a > the symplectic leaf M/3_q, 
is a submanifold of M of dimension 2(A^— 1), there are at most A^— 1 independent 
integrals in involution on Mp.a- For any given (6, a) £ M/3,q let Trp,a denote the 
orthogonal projection Tf,^aM Tb^a^i3,a- Then the gradient of the restriction 
/n|]vt3,c of In to M/3,Q {I < n < N—1) IS given by 1^(3,0.^ b,aln- By Lemma Fr2l the 
vectors [iip^a^h.alnjneK are linearly independent. As M/3_q \ ^k=i^k is dense 
in JAp^a, it then follows that (/n|3vt^j „)i<n<A^-i ^'''^ functionally independent. 
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Finally, as Ci and C2 are Casimir functions of {•,•},/, it follows that for any 
(6, a) e Mp^a 

= {^bMln, J'^b,alm) — {In,l7n}j = 0, 

i.e. the restrictions In\Mf, „ of /„ (1 < n < — 1) are in involution. □ 

A Proof of Lemma 13.71 

In this Appendix, we prove Lemma [3.71 It turns out that the proof in ([1], 
p. 601-602) of the special case where the parameter a in ([1]) equals 1 can be 
adapted for arbitrary values. 

Proof of Lemma \3.7\ Let (6, a) be an arbitrary element in M and I < n < N-^l. 
First note that /„ ~ ^ /a2^"^^ arcoshj iA(A)| dA and use ^arcosh(t) = ^ 
to obtain 

^ ^A2„ + i ^|A(A)|/2 ^ 

- / , dt dX. 



Since the integrand of the inner integral is nonincreasing, we estimate it from 

2 



below by its value at This leads to 



/„ > i dX. (117) 



TT 7,,„ ^A(A)| 

We will show below that for X2n < < ^2n+i 



V|A(A)| - 2 ^ VX^X^^X2n+i - X 
V|A(A)| +2 " A2JV - Ai 

We then substitute (IllSp into the integral pi7p and split the integration interval 
into two equal parts. 



, ^ 2 1 Z"^" \/A — A2n V A2n+1 " A 

^"-~~\ r ^ — \ 

71" A2N — Ai Jx^^ Mn — Ai 

where t„ = (A2„ + A2„+i)/2. For A2n < A < we estimate the quantity 
A2n+i — A from below by 7n/2, yielding 

/„ > ^ . r a/¥v/a^ dx = ^ 



TT A27V — Al 7;^^^ V 2 ^ ^" 37r(A2Ar — Al) ^"^ 

It remains to verify (|118p . Recall that A2n and A2„+i are either both periodic 
or both antiperiodic eigenvalues of i. If A2ti and A2n+i are periodic eigenvalues, 
we have A(A) > 2 for A2„ < A < A2r!+i, i.e. |A(A)| = A(A). In order to make 
writing easier, let us assume that N is even - the case where N is odd is treated 
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in the same way. Then by Ai and A2Ar are periodic eigenvalues of L and 
thus for any < A < A2n+i, the left side of (|118p can be estimated from 
below by 



VA(A) - 2 ^ / (A-A2„)(A-A2„+i) ^ ^ VA - A2„v/A2n+i - A ^ 

V'A(A) + 2 V (^"^2)(A-A2W-l) ' " A2Ar-Ai ' ' 



where 



R^RW = J^---^-^Y^---^^^- (119) 

y A — A3 A — A2n-1 A2ri+2 ^ A A2JV-2 ^ A 

As each of the the fractions under the square root in (I119P can be estimated from 
below by 1, for any A2ri < A < A2n+i it follows that i?(A) > 1 on [A2„,A2ri+i], 
leading to the claimed estimate (IllSp . □ 

B Proof of Theorem Q 

In this Appendix we prove Theorem l3 . 91 using estimates derived in J]. Let (6, a) 
be in M/j.q with /? S M and a > arbitrary. 
To show Theorem 13.91 we need the following 

Proposition B.l. For any {b,a) G M-fj^a with /3 e M, a > arbitrary and any 
l<n< N, 

A2„(&,a)-A2„_i(6,a) < — . (120) 



Before proving Proposition IB . II we show how to use it to prove Theorem l3.9l 
Proof of Theorem \3.9l We begin by adding up the inequalities ([60|) and get 

N-l /N-1 \ 

^7,'.<37r(A2Ar-Ai)^/„ . (121) 

Note that 

N-l N 
^2N — Ai = ^ 7„ + ^(A2„ — A2„-l). 
n— 1 n— 1 

By the estimate of Proposition IB . II we get for any (6, a) G Mp^a 

N~l 

\2N - Ai < 27ra + ^ 7„ 

n=l 

which we substitute into (I12ip to yield 

N-l /N-l \ /N~l \ /N~l \ 

5:7^<6a.^K:/„ +3. K:7jK:/j. 

n=l \n=l / \n=l / \n=l / 
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Using the inequality 

2ab< e'^a'^ + ^b'^ (a,6eR,e>0) 

with a = Y.n=i ln,b^ Y.n=i ^n, and = 3.^(^n-i) ^ ^"^^ S^*^ 

N-l /N-1 \ o / 1 /AT-l ^ ^ 



(l]n=i^ 7«) < - 1) (Z]n=i^ 7n) ) ouc thcn concludcs that 

Af-l /N-l \ 2 /^-l \ ^ 

iE^^^STT^aK]/,, +^(iV-l)K:/J , (122) 



n— 1 \n— 1 / \n—l / 

which is the claimed estimate (IFTI). □ 



To prove Proposition IB. II we first need to make some preparations. Note 
that for an element of the form (6, a) — (/SIat, a^N) one has, by Lemma [ 



A2n(/3lAf, alAr) - A2„-i(/3lAr, alw) = 2a|^cos 



(n — 1)tt mr\ 
cos — I 



(2n — l)7r TT 

= Aa sin — — sin — — 

2N 2N 

2'Ka 

Hence to prove Proposition IB . 1 1 it suffices to show that for any (5, a) e M/j^q, 
and any 1 < n < 

A2n(fe, a) - A2„-i(&, a) < A2„(-/9lAr, alAr) - A2„-i(-/31jv,alAr)- (123) 
To this end, following [Kj (cf. also [7j), we introduce the conformal map 

6{X):=[-ir f ^P^ d,, (124) 



where the sign of the square root is chosen such that for /i < Ai , ■\/4 — A^(/i) has 
positive imaginary part. It is defined on the upper half plane U ;= {Imz > 0} 
and its image is the spike domain 



n(&, a) -.^ {x + iy: Q < X < Ntt, y > Q} \ |J r„ 

11=1 

where for 1 < n < iV — 1 , T„ denotes the spike 

T^:=\nn + it:Q<t< arcosh ( (_l)JV+"^^M j 
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To see that 5{U) = n{b,a), note that for any (6, a) G M and X E U the 
discriminant A(A) and the function S{X) are related by the formula 

A(A) 2(-l)^ cos 5(A). (125) 

To prove (|125l) . recall that for — 1 < < < 1, one has ^arccosi = -^j=^. 
This formula remains valid for any t in C \ ((— oo,— 1] U [1,cxd)). Thus 



5{\) = (-1) 



- t^M^d,= r:^arccosf(-l)-^) 



Since by dH]), A(Ai) 2(-l)^, we then get 

5{\) = arccos f(-l)^^^) ^ = arccos ( (-1)"^^^ ) , (126) 



leading to formula (|125[) and the claimed statement that 5{U) = n{b,a). 

The map S can be extended continuously to the closure {Imz > 0} of the 
upper half plane. This extension, again denoted by S, is 2-1 over each nontrivial 
spike r„ and 1-1 otherwise. Since the n-th spike T„ is the image under 6 of 
the n-th gap (A2„,A2„+i), all spikes are empty iff all gaps are collapsed. By 
Lemma [2.61 all gaps are collapsed for {b,a) = (— /31jv, alAr), hence f2(6, a) C 
n(— /31jv, alff) for any {b, a) G Mp^a- 

Note that 

A2n(6,a)-A2„-i(6,a)=ri(n7r-)-ri((n-l)7r+)= / u^"\S{X))dX, 



(127) 

where u'") : n{b, a)^R, u'") {z; b, a) is the harmonic measure of the open 
subset {{n — l)7r, nn) of dil{b, a) (see e.g. [Bl for the notion of the harmonic 
measure). 

We need two lemmas from complex and harmonic analysis, respectively. 

Lemma B.2. For (6, a) = (— /31jv, aljv) with /3 e M, a > arbitrary, the map 
(5(A) defined by \124^ is given by 

S{X)^N arccos ■ (128) 

For arbitrary (5, a) in M^,a and ^ G M, the following asymptotic estimate holds 
as rj oo 

5(e + "/) -iV arccos +0{r]-^), (129) 

locally uniformly in ^. 

Proof of Lemma \B.SX In view of the formulas ((5T|) and ((52|) for the fundamen- 
tal solutions yi and y2 for (&, a) = {— [31 m , Oil n) ^ the discriminant A(A) = 
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A(A, — /31jv, alw) is given by 

A(A) = yi{N,X)+y2{N + l,X) 

sin(p(A^-l)) sin(p(iV + l)) 
sin p sin p 

= 2cos{pN), 

where tt < p < 27r is determined by cos p ~ . Hence 



A(A) = 2T^. ( ^) (130) 



2a 

where for any z £ U, 

Tjy (z) = cos{N aiccosz). (131) 

Actually, Tm{z) is a polynomial in z of degree iV, referred to as Chebychev 
polynomial of the first kind. Substituting (|130p into ()126p . we obtain 



Af^] ^ arccos f (-l)^r^ 



(5(A) = arccos l^(-l) ^ ^ ^ 

The claimed identity (|128p now follows from the elementary symmetry 

Tw(z) = (-l)^TAr(-z) VzeC. 

Now let (6, a) G M^,a- The asymptotic estimate (I129P follows by compar- 
ing the polynomials A (A) correponding to (6, a) and the one corresponding to 
(— /JIat, aljv)- By (fT7|) (and the discussion following it), in both cases. 



A(A) = a-^A^ (1 + 7V/3A-1 + 0(A-2)) as |A| 
This implies that 

A,,,(A) = A_^i„,„i„(A) • (1 + 0(A-2)), 
hence by and ([^0]) 



oo. 



5,,,(A) = arccos (^(-1)^^^^) 



= arccos f A_^i„,„i„ (A) • (1 + 0{X-^)) 

= arccos (^(-l)^rAr (^^) • (1 + 0(A-^)) 

= arccos (^Tat ' (1 + 0(A-2)) 

Substituting formula (|13ip for T/v, one then concludes that 

6b^a{X) = arccos ^cos arccos ' + '^('^~^)) 

= A arccos ('-^^'j +0(A-2), (132) 
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where in the last step we used that arccos z = —i log(z + i \/l — z"^) for any z in 
C\ ((-oo,-l] U [l,oo)). □ 

Lemma B.3. Let u : 51(6, a) ~^ R be a bounded harmonic function such that the 
nontangential limit of u{z) on dfl{b, a) has compact support, and let U{X) := 
u{S{X)), where (5(A) is the function defined by \124^ . Then for almost every 
i e M, the limit U{t) := lim.,,^o U{t + irj) exists and is integrable, and 

f°° / x\ f Nt: \ 

/ U{t)dt = Km 27ra sinh — it hix] . (133) 

J-oo ^NJ \ 2 ) 

Proof of Lemma \B.3[ Again by Fatou's theorem, for a.e. t the (nontangen- 
tial) limit lim^^o U{t + irf} exists, since C/ is a bounded harmonic function on 
{Im {z) > 0}. Since u is bounded on n(b,a) and its nontangential limit to dil 
has compact support, U{t) is bounded and compactly supported and thus in 
particular integrable. For X ^ + irj, one then has the Poisson representation 

and by dominated convergence we conclude that 

U{t)dt^ lim TnjU{^ + ir]). (135) 

(In particular, the limit in the latter expression exists.) In order to compute 
the right hand side of (|135p , let ^ + irj be given by 



£, + iri = I — — h ix j 



for X sufficiently large. Then J7(^ + ii]) = u + i.x), and by ()132p . it follows 
that 



TT i f iC + f3) + iv 

— X = arccos ' 

2 N 



^JA+fl±2!lyO{{^ + ^rJr') (x^oo). (136) 



Taking the cosine of both sides of (|136p . multiplying by —2a and using that 
cos(|- + it) = —i sinhi for t G K, we obtain 

2iasinh^ = {{^ + 13) + i?]) {1 + O {{£. + iTj)-^)) (a; -> oo). 



Hence, as a; ^ oo, 



e = 0(l), 77 = 2a sinh ^ + 0(1). (137) 



Substituting p37p into (|135p leads to the claimed formula p33p . □ 
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Proof of Provosition \B.l\ Let (6, a) G M/j.q. Besides the harmonic measure m^") 
of the set E := {{n — 1)7t, mr) C dfl{b, a) we also consider the harmonic measure 

u^p l, of E C di}{~l3lN,alN)', note that 

Q{~(3lN,alN) ^ {x + iy\0 < X < N-k, y > 0} 

and hence i}{b, a) C ri(— /Jl^r, al^v)- According to |6j, both m^") and satisfy 
the hypotheses of Lemma IB. 31 Let us recall the monotonicity property of the 
harmonic measures u{z,E,n) with respect to (see e.g. [6]): If ili C 
E C d^li n d^l2, and u{z,E,ni) [i = 1,2) denotes the harmonic measure of 
E C 9f2i, then for any z £ fti, u{z, E,ni) < u{z, E,rt2)- Apply this general 
principle to fli :— Q{b,a) and Q2 '■= ^{—P^n,0£1n) to get 

< (138) 

Writing t/(")(A) := u^"\S{X)) as weU as C/^"](A) := and combin- 

ing dnH), dnSl), and (UnSl), we conclude that 



X2n{b,a)~X2n-i{b,a) ^ U^"\X)dX 

= lim 27rasinh I — I u'"^ 



/Nn 
/Nn 



< lim 27rasinh( — | u''"'' 



= A2n(-/3l7V, alw) - A2„-l(-/3lAr, al^v)- 

This completes the proof of the estimate (I120p and therefore of Proposition 

[bH □ 
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